Analysis of two-dimensional high-energy photoelectron momentum distributions in 
single ionization of atoms by intense laser pulses 
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We analyzed the two-dimensional (2D) electron momentum distributions of high-energy photo- 
electrons of atoms in an intense laser field using the second-order strong field approximation (SFA2). 
The SFA2 accounts for the rescattering of the returning electron with the target ion to first order 
and its validity is established by comparing with results obtained by solving the time-dependent 
Schrodinger equation (TDSE) for short pulses. By analyzing the SFA2 theory, we confirmed that 
the yield along the back rescattered ridge (BRR) in the 2D momentum spectra can be interpreted as 
due to the elastic scattering in the backward directions by the returning electron wave packet. The 
characteristics of the extracted electron wave packets for different laser parameters are analyzed, 
including their dependence on the laser intensity and pulse duration. For long pulses we also studied 
the wave packets from the first and the later returns. 

PACS numbers; 32.80.Rm, 32.80.Fb, 42.50.Hz 
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INTRODUCTION 

Much of our knowledge of strong field, multiphoton 
ionization processes conies from the study of above- 
threshold ionization (ATI). In this process an atom ab- 
sorbs more photons than the minimum number required 
for the ejection of an electron. By measuring the en- 
ergy of the photoelectrons, the characteristic ATI spec- 
tra are peaks separated by the photon energy, with the 
peak positions shifted by the ponderomotive potential. 
For long pulses, additional substructures due to Freeman 
resonances [1] appear. For shorter pulses, of durations 
of the order of 10-20 fs or less, calculations have shown 
that the major ATI peaks are accompanied by subpeaks 
which have been attributed to the rapidly changing pon- 
deromotive potential ^2, 3]. 

Experimentally more information on ATI electrons can 
be determined by measuring the angular distributions or 
the 2D momentum distributions. For example, the na- 
ture of Freeman resonances are associated with individual 
Rydberg states by the number of lobes in their angular 
distributions 0] ■ For short pulses in the tunneling ioniza- 
tion regime, recent experiments have shown ubiqui- 
tous fan-like structures in the 2D momentum spectra for 
low-energy electrons. These fan-like structures are now 
also well-understood and they are due to the long-range 
Coulomb interaction between the tunnel ionized electron 
and the target ion it has left behind 0, Q • These fan-like 
structures do survive the integration over the laser focus 
volume, as shown in recent comparison between theoret- 
ical calculations with experiments Q . 

High-energy photoelectrons have been measured previ- 
ously for atomic targets, both in energy and angular dis- 
tributions, using longer pulses at lower intensities. The 



energy spectra above Wp, where Up is the ponderomo- 
tive energy, have been found to vary rapidly with small 
changes in laser intensities |l3, 11 1 when laser pulse du- 



rations are larger than 20 to 30 fs. Different theoretical 
models have been used to interpret these phenomena, but 
the conclusions are still tentative so far 12, 13, 14, T^. 
While single ionization of an atom in an intense laser 
field has been calculated accurately by solving the time- 
dependent Schrodinger equation (TDSE) directly in the 
past two decades, the photoelectrons in the high-energy 
region, especially their momentum distributions, have 
not been investigated. This is not surprising since the 
photoelectron yield drops rapidly with the electron's en- 
ergy. On the other hand, it is known that rescatter- 
ing plays a major role in many laser-atom interaction 
phenomena, including the generation of high-energy ATI 
electrons. Understanding the nature of these high-energy 
photoelectrons, especially their 2D momentum distribu- 
tions, may help to shed new light on the rescattering 
process itself. 

Recently we have initiated a careful study on the 2D 
momentum distributions of high-energy photoelectrons 
by directly solving the TDSE of a one-electron atom in 
a short laser pulse. By focusing on photoelectrons that 
have been backscattered we have been able to extract 
the elastic scattering cross sections by free electrons from 
the calculated laser-generated photoelectron spectra. For 
molecular targets this has the important significance that 
it offers the possibility of using infrared lasers for dy- 
namic chemical irnaging with temporal resolution of a 
few femtoseconds [l6[. For short laser pulses, by study- 
ing the dependence of the 2D momentum spectra on the 
carrier-envelope phase, it also offers the possibility of di- 
rectly characterizing the electric fields of few-cycle pulses. 
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To obtain 2D momentum electron spectra accurately 
in the high-energy region, say, up to about 10 Up, cal- 
culations have to be done very carefully. For pulses as 
short as a few femtoseconds and intensities with Keldysh 
parameters 7 17| close to one, we have been able to ob- 
tain accurate 2D momentum spectra from solving TDSE 
directly. Here 7 — \Jl-pl (2C/p), with /p the binding en- 
ergy of the electron, and C/p = /nia2;/(4ci;^), where l^aax 
is the peak intensity and lo is the carrier frequency, of 
the laser, respectively. (Atomic units are used through- 
out unless indicated otherwise.) To extend such TDSE 
calculations to high laser intensities or long durations is 
computationally more challenging. We thus seek to ex- 
amine the predictions based on the second-order strong- 
field approximation (SFA2) [H, [H, [i^, [2l|, [H, where 
the rescattering of the returning electron with the target 
ion is included. To establish the validity of SFA2, we first 
compare its results with those from solving the TDSE 
for short pulses. The SFA2 is then used to analyze high- 
energy photoelectron momentum distributions for vari- 
ous laser parameters. From such analysis, we show that 
electrons on the so-called back rescattered ridge (BRR) 
can be identified. The electron yield along the BRR can 
be used to extract elastic scattering cross sections of free 
electrons by target ions, as well as the wave packet of 
the returning electrons. Using SFA2 allows us to analyze 
the electron wave packets for longer pulses and higher 
intensities. 

In the next Section we first summarize the first- and 
second-order amplitudes of the strong field approxima- 
tion. The results from using this theory for the predic- 
tion of the 2D momentum spectra are shown in Section 
III. The conclusion is given in Section IV. 



THEORETICAL MODELS 

To study atoms in an intense laser field, there are only 
two general theoretical tools at present. One is to numer- 
ically integrate the TDSE and the other is based on the 
SFA. Both approaches have their limitations. For TDSE, 
to obtain high precision in small quantities such as the 
momentum distributions of high-energy electrons poses a 
computational challenge. This challenge is less severe for 
short pulses where the electron is confined to a smaller 
box such that accurate solution of the TDSE is possible. 
For longer pulses and higher intensities, high accuracy 
becomes harder to achieve. One has to check carefully 
the possible effect of reflection from the box boundaries 
as well as the convergence of the basis set used. While 
reflection can be avoided by matching the solution in- 
side the box to the outside region where it is expanded 
in terms of Volkov states ji^l , such an approach requires 
large basis set in both the inside and the outside regions, 
and thus the method is also very computationally inten- 
sive. Thus in the TDSE calculations we limit ourselves 



only to short pulses. The numerical method we used for 
solving TDSE has been presented in our previous works 
0, In this paper we only present the details of the 
SFA2 used here. 

An exact expression for the probability amplitude of 
detecting an ATI electron with momentum p can be writ- 
ten formally as 



/(P) 



lim 



dt'(vl/p(i)|[/(i,t')^.(i')l*o(i'))-(l) 



Here U(t, t') is the time-evolution operator of the com- 
plete Hamiltonian 



H{t)=Ha+H,{t) 



where 



Ha 



y(r) 



is the atomic Hamiltonian, and 



(2) 



(3) 



(4) 



is the laser-electron interaction in the length gauge and 
the dipole approximation. The linearly polarized electric 
field E(i) of the laser pulse along the z axis is given by 



E(t) = Eoa{t) cos{ujt + (j))z 



(5) 



where (j) is the carrier-envelope phase. The envelope func- 
tion a{t) is chosen to be 



ait) 



cos 



(6) 



for the time interval (— T/2,T/2), and zero elsewhere. 
In this paper, T is defined as the (full) duration of the 
laser pulse which is 2.75 times of the FWHM (full width 
at half maximum) and the carrier-envelope phase </> is 
set as zero. The functions vE'p(t) and \l/o(i) are scatter- 
ing state with asymptotic momentum p and the ground 
state, respectively, of the atomic Hamiltonian Ha- The 
time-evolution operator U (t, t') satisfies the Dyson equa- 
tion 

U{t,t') = UFit,t')-i f dt"UFit,t")VU{t",t') (7) 



where UF{t,t') is the time-evolution operator for the 
Hamiltonian of a free electron in the laser field, which 
is 



1. 



r -E. 



(8) 



The eigenstates of Hp (t) are the Volkov states 

|xp(i)) = |p + A(i))expH^pW] (9) 
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with the action 



Spit) ^ I f_Jt' [v + A'it)]' 



(10) 



The vector potential of the laser field E(<) is denoted by 
A{t), and |k) is a plane wave state 



(rlk) = 



1 



■ exp (ik • r) . 



(27r)3/2 

The Volkov time-evolution operator is 

UAt,t')^ f dk|xk(t))(xk(0 



(11) 



(12) 



By approximating U{t",t') on the right-hand side of 
Eq. (7) by UFit",t'), and {^pit)\ in Eq. © by (xk(i)|, 
the ionization amplitude may be expressed as 



f(2) 



(13) 



where the first term 



dt(xp(t)|ff.(t)|*o(t)) (14) 

corresponds to the standard SFA. This term will be called 
SFAl for the present purpose. The second term is the 
SFA2, 

/OO nt n 

dt dt' dk{xp{t)\V\x],{t)) 
-oc J — OO J 



x(xk(t')|i^.(i')l*o(t')) 



(15) 



which accounts for the first-order correction by the 
atomic potential. This expression can be easily under- 
stood by reading it from the right side. The electron is 
first ionized at time t' by the laser field, ft then propa- 
gates in the laser field from t' to t where it is rescattered 
by the atomic potential V into a state with momentum 
p. Note that in this approximation, the interaction of 
the electron with the atomic potential is treated up to 
the first order only. 

The evaluation of the matrix elements are illustrated 
below for the hydrogen-like atoms where the ground state 
wavefunction takes the form 



*is(r) = 2^3/2 cxp(-Zr)roo(f) 



(16) 



where Z is the charge of nucleus. The rescattering po- 
tential V{v) is a pure Coulomb potential 



V(r) 



(17) 



In the numerical calculations, the Coulomb potential is 
replaced by Yukawa potential with a damping parameter 
a to avoid the singularity in the integrand in f^'^\ 



V{v) = V{v)e- 



(18) 



This introduces a very weak dependence of the magni- 
tude of /^2) on the value of a, but not the shape. Here, 
we chose a — 1.0. 

After performing integration over space coordinates 
analytically, the amplitudes /'■^•' and /f^^) become 



^(1) ^ _^^5/2 /■ dtexp[z^p(i)]exp(i/pi) 

^ J —OO 



and 



{Z2 + [p + A(t)P} 
AV2 



(19) 



dt 



dt' / dk 



exp {-z[S'k(t) - Spit)]) exp[iS'k(t')] exp(i/j,i') 
1 E(i') • [k A(t')] 



(k-p)2{^2 + [k + A(t')]2}' 



(20) 



respectively. The evaluation of the first-order amplitude 
(fTO)) is straightforward. The second-order amplitude (^0]) 
consists of fivefold integration. We used saddle point 
approximation for the integration with respect to k, as 
proposed in Lewenstein et al l^l ■, to reduce it to a twofold 
one. The saddle point is calculated with respect to the 
quasiclassical action only, so that 



ks(i,i') = - 



1 



t-t' 



dt"A{t" 



(21) 



The result of saddle point integration for the fivefold in- 
tegral ((20|) is obtained by setting 



k = ks(i,t'), 



and substituting 



dk 



2lT 



e + i{t~t') 



3/2 



(22) 



(23) 



for the integration over k. Here, e is an arbitrary small 
parameter introduced to smooth out the singularity in 
(1221) ■ Its value is taken to be of the order of 0.1. 

The momentum distribution of the emission of an elec- 
tron of energy E — in the direction of p is given by 



dEdp 



(24) 



The form of / in allows us to identify the contribu- 
tion from each individual terms SFAl and SFA2, respec- 
tively. 

For a linearly polarized laser field, the system has cylin- 
drical symmetry for the cases considered here. As a re- 
sult, the two-dimensional momentum distribution is de- 
fined by 



d^P 



dEde dEdp 



2'Kp sin t 



(25) 
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where integration over the azimuthal angle has been 
carried out and 6 is the angle between the polarization 
axis of the laser field and the direction of the ejected 
photoelectron. By integrating over 6 in equation (|25p . 
we obtain the energy spectra 



dP 
dE 



dEdO 



de. 



(26) 



RESULTS AND DISCUSSION 
Validity of the SFA2 

First we establish the validity and the limitation of the 
SFA2. In Fig. 1(a) we show the total ionization prob- 
ability vs electron energy for a hydrogen atom ionized 
by a 5-cycle (full duration) laser pulse with the wave- 
length of 800 nm and at the peak intensity of 1.0 x 10^'' 
W/cm^. For this case, the ponderomotive energy is 6 
eV. Within the perturbation approach we note that the 
yield is dominated by SFAl for low energy electrons. For 
energies greater than about 4J7p, SFA2 becomes domi- 
nant. Interference between SFAl and SFA2 is important 
only in a small energy region. Note that the actual prob- 
ability obtained from SFAl severely underestimates the 
total ionization yield as obtained by the TDSE. This un- 
derestimate is present also in SFA2 since the same ma- 
trix element for the initial ionization of the atom is also 
used in SFA2. Here we are interested in the high-energy 
region, we thus renormalize the electron energy spectra 
to the TDSE results at higher energies. Note that after 
normalization, see Fig. 1(b), the electron energy spectra 
from the two theories are in good agreement for energies 
above blip. This comparison also proves that rescattering 
plays a major role in the high-energy ATI electrons. 

We next consider the 2D electron momentum spectra. 
Previously theoretical calculations tend to focus on the 
angular distributions at specific ang les, especially along 
the laser polarization direction 25|, |26|. Our goal, in- 
stead, is to examine the global 2D momentum spectra 
in the high-energy region and whether the spectra can 
be described by the SFA2. In Fig. 2 the 2D momentum 
spectra calculated from SFA2 and from TDSE are shown. 
The horizontal and the vertical axes are the electron mo- 
menta parallel and perpendicular to the laser polariza- 
tion, respectively, on a plane containing the polarization 
axis. Since the electron yield drops very rapidly with 
energy (see Fig. 1), to make the high energy part visi- 
ble, the 2D momentum distributions in each frame have 
been renormalized such that the total ionization yield at 
each electron energy is the same. At first glance, clearly 
the two frames look very similar in the high-momentum 
region. (The difference in the small momentum part is 
not important since SFAl is dominant there.) Similar 
agreement has also been observed for lasers of different 
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FIG. 1: (Color online) Electron energy spectra for atomic hy- 
drogen by a 5-cycle laser pulse with the wavelength of 800 nm 
at the peak intensity of 1.0 x lO" W/cm^(a) From TDSE, 
first- (SFAl) and 2nd-order (SFA2) theory and the coherent 
sum of SFAl and SFA2; (b) the sum of SFAl and SFA2 but 
renormalized to the high-energy part of TDSE. Electron en- 
ergy is in units of Up, the ponderomotive energy. 



intensities and wavelengths. 

The striking features of Fig. 2 are the two half cir- 
cles on each side of the origin. The center of each cir- 
cle is shifted along the p\\[= Pz) axis. We call these 
circular rings back rescattered ridge (BRR), represent- 
ing electrons that have been rescattered into the back- 
ward directions by the target ion. Based on the results 
from TDSE calculations, we have confirmed [l^ that the 
electron yields on the BRR can be interpreted as elas- 
tic scattering of the returning electrons by the target ion 
potential. Since the main features of the BRR are also 
reproduced in the SFA2 we seek to provide theoretical 
grounds for this interpretation. 

Evidence of these BRR electrons had been seen ob- 
served previously. Using 50-ps, 1.05-Atm linearly polar- 
ized laser pulses, Yang et al. [27[ observed unexpected 



narrow lobes in the angular distributions at approxi- 
mately 45° off the polarization axis in the high order 
ATI spectrum around 9Up in single ionization of Xe and 
Kr atoms. The origin of the sidelobes has been analyzed 
by Paulus et al. [28| in a two-step classical model and 
by Lewenstein et al. [l^l using a quasiclassical analysis. 
Similar theoretical analysis of these narrow sidelobes has 
been made by Dionissopoulou aL [2^. However, direct 
connection of these narrow sidelobes with the backscat- 
tering of the returning electrons has not been established 
at the quantitative level so far. 



Analysis of the high-energy 2D momentum spectra 

To identify that rescattering is responsible for the high- 
energy BRR electrons, in Fig. 3 we show the renormal- 
ized 2D momentum spectra for 2-cycle and 3-cycle 800 
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FIG. 2: (Color online) Photoelectron 2D momentum distri- 
butions parallel (py = p^) and perpendicular (py) to the laser 
polarization direction for atomic hydrogen by a 5-cycle laser 
pulse with the wavelength of 800 nm at the peak intensity of 
1.0 X 10" W/cm^ (a) SFA2 and (b) TDSE. 

nm pulses at peak intensity of l.Ox 10^^ W/cm^, with the 
corresponding electric fields and vector potentials plotted 
on the right side. The data in the 2D momentum spec- 
tra are integrated over the azimuthal angle around the 
polarization vector, as in (j25p . This presentation forces 
the distribution to go to zero on the p|| axis. 

Consider the 2-cycle pulse, from Fig. 3(b) the elec- 
tric field E{t) reaches local maximum at time t = — 0.4r, 
0, and 0.4r, where r = 2tt/uj is the period of the laser. 
Based on classical theory, electrons that are " born" near 
the peak of the electric field (more precisely, at the phase 
angle of 17° after the peak) will be driven back to the tar- 
get ion about three quarters of a cycle later. Consider an 
electron that was born at the time tb near — 0.4t where 
the electric field strength is maximum. We will choose 
the convention that positive z is the right side and nega- 
tive z the left side. The released electron is in a negative 
electric field so it will first moves to the right side. After 
the electric field changes to positive, the electron will be 
decelerated and may be driven back toward the target. 
Simple classical calculation shows that it will return to 
the origin at time U near 0.25r where the electric field 
returns to zero. If this electron is scattered in the forward 
direction, it will be decelerated by the subsequent nega- 
tive electric field for t > tr and ending up with low energy. 
On the other hand, if the electron is backscattered elasti- 
cally, it will be further accelerated by the laser's electric 
field and ending up with high energy. This effect on the 
electron momentum is obtained by adding the instanta- 
neous vector potential = A.{tr) to the canonical mo- 
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FIG. 3: (Color online) Photoelectron 2D momentum distribu- 
tions parallel (py = Pz) and perpendicular {p± — -y/pi + Py) 
to the laser polarization direction for atomic hydrogen by laser 
pulses at the peak intensity of 1.0 x 10^* W/cm^ with the 
wavelength of 800 nm for durations of 2 cycles and 3 cycles, 
respectively. The corresponding laser fields and vector poten- 
tials are shown on the right for the analysis of the rescattering 
mechanism. 



mentum p, analogous to the "streaking" of an electron 
generated by an X-ray attosecond pulse in the presence 
of a femtosecond IR pulse [2^ . The shift of the center of 
the BRR is a measure of the magnitude of A^ , and the 
radius of the circle is related to the maximum electron's 
returning energy 3.17Up. In Fig. 3(b), electrons that are 
born near t — 0, would return at about t — 0.75t where 
the vector potential and the electric field are near zero, 
thus no high-energy electrons would emerge. This ex- 
plains why there is no BRR on the left side of Fig. 3(a) 
for the 2-cycle pulse. 

Consider next the 3-cycle pulse, the right BRR is due 
to ionization occurs near tb and returns near tr as before. 
The BRR on the left is due to ionization near t'f^ and 
returns near t'^.. Since the vector potential at t'^. is smaller 
than that at tr, the shift of the center and the radius of 
the circle of the left BRR are smaller. 

One advantage of the SFA2 is that it allows us to iden- 
tify the born time and the returning time directly. To 
analyze the right-side BRR in Fig. 3(c), we set "win- 
dow functions" such that the born time is restricted to 
[— 0.5r, — 0.4r] for t' and the return time interval [0, 0.75r] 
for t in (|20)) . as shown in Fig. 4(b). The resulting 
2D spectra shown in Fig. 4(a) is similar to the right 
side of the 2D spectra in Fig. 3(c), confirming that the 
right BRR is from electrons born near tb- Similarly, in 
Fig. 4(c), the 2D momentum distributions are presented 
to identify the left-side ridge by setting the born time 
interval [0,0. Ir] for t' and the returning time interval 
[0.5r, 1.5r] for t in as shown in Fig. 4(d). The re- 
sulting left BRR is also similar to the one in Fig. 3(c). 
These confirm our interpretation of the origin of each 
BRR, in terms of its born time and returning time. 



FIG. 4: (Color online) Analysis of 2D momentum distribu- 
tions as in Fig. 3 for the 3-cycle pulse, for the identification of 
the tunneling time and rescattering time, for the right- and 
left-side BRR electrons. See text. 



The position of the BRR in the 2D momentum space 
can be expressed as 

p = -A, + (27) 

where p^ is momentum vector measured from the center 
of the circle which is located at — and hence pr is the 
radius of the BRR. The projection of the photoelectron 
momentum in the parallel and perpendicular directions 
are 



P|| = PC0S9 = —Ar — PrCOsOr 

p± = p sin 9 = Pr sin Or 



(28) 



where 9r is the backscattering angle, ranging from 90° to 
180°. 

Actually, the BRR in the 2D momentum distribution 
is formed by the cutoff peak in the angle-resolved energy 
spectra. To check the accuracy of the interpretation rep- 
resented by (P7)) . we "measured" the center position and 
the radius of the right-side BRR as \Ar\ = 0.93, Pr = 1.18 
in the 2D momentum distributions shown in Fig 2. We 
chose 9r = 180°, 160°, 140°, and 120°, which correspond 
to 6* = 0°, 11.2°, 22.5°, and 34°, with the corresponding 
energies of 2.19, 2.13, 1.94, and 1.65, respectively. In Fig. 

5, we show the angle-resolved energy spectra at these an- 
gles with the predicted electron energies by arrows. The 
arrows indeed are located at the peak positions of the 
angle-resolved energy spectra. 

According to the previous paragraph, the BRR de- 
pends on the vector potential of the pulse only. In Fig. 

6, we show the renormalized 2D momentum spectra for 
single ionization of hydrogen atom by 5-cycle laser pulses 
with fixed Keldysh parameter 7 — 1.07. The wavelengths 
are 400, 600, 800, and 1000 nm, and the corresponding 
peak intensities are 4.0 x lO^'', 1.78 x lO", 1.0 x lO^'', and 
6.4 X lO^'^ W/cm^, respectively. It can be seen from Fig. 
6 that the BRR's exist for all the cases and the 2D spec- 
tra appear to be very similar to each other. It can also 
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FIG. 5: (Color online) Angle-resolved energy spectra for 
atomic hydrogen by a 5-cycle laser pulse with the wavelength 
of 800 nm at the peak intensity of 1.0 x 10^'' W/cm^, for four 
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FIG. 6: (Color online) 2D momentum distributions for atomic 
hydrogen by 5-cycle laser pulses with the same Keldysh pa- 
rameter 7 = 1.07, but for wavelengths of 400, 600, 800, and 
1000 nm, respectively. 



be seen in Fig. 6 that the BRR breaks into sub-rings and 
the width of the BRR becomes narrower with increase of 
wavelength. These sub-rings are not investigated here. 



Returning electron wave packets 

By identifying the BRR electrons as due to the 
backscattering of the returning electrons by the target 
ion, the yield along the BRR should then be proportional 
to the differential elastic cross sections of electrons by 
the target ion. Within the SFA2, the elastic scattering is 
treated in the Born approximation, where the scattering 
amplitude is proportional to the Fourier transform of the 
ion potential 



V{q) = J exp(iq • r)\/(r)dr, 



(29) 
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FIG. 7: (Color online) Right-side BRR for atomic hydrogen 
by a 5-cycle laser pulse at the peak intensity of 1.0 x 10^'* 
W/cm^ with the wavelength of 800 nm. (a) Comparison of 
the first Born elastic scattering cross sections with the angular 
distributions on BRR calculated from SFA2; (b) Comparison 
of the electron wave packet from SFA2 with that from TDSE 
at 6r = 160° ; (c) The wave packets of the returning electrons 
are shown to be independent of the scattering angles. 



here q is the momentum transfer which is related to the 
rescattering angle 9r and the radius pr of the BRR by 



2pr sm{9r/2). 



(30) 



Fig. 7(a) shows the comparison of the first Born elec- 
tron elastic scattering cross sections with the angular dis- 
tributions on the right-side BRR in the 2D momentum 
spectra shown in Fig. 2(a). For the comparison pre- 
sented in Fig. 7(a), the radius of BRR is taken to be 
Pr = 1.2. It can be seen that for 9r > 100°, the differen- 
tial elastic scattering cross sections are very close to the 
photoelectron angular distributions along BRR, indicat- 
ing that elastic scattering of electrons by the parent ion 
can be extracted from the angular distribution along the 
BRR. 

Expanding the rescattering idea further, we then ask 
if it is possible to treat the electron yield on the BRR 
as due to the backscattering of a returning electron wave 
packet. To test this idea, we write 



dEdp 



I c^— cons 



HPr,Or)\V{q)\' 



(31) 



On the left side in |3T|) . the angular distribution is ob- 
tained on any plane containing the z axis in the momen- 
tum space due to the cylindrical symmetry. Note that 
|y(q)P is the elastic differential cross section for each 
ion by an electron with energy Er = Pr/2 (see (15(1)) ) 
while the actual energy of electron in the laser field is 
E = [—Ar + Pr]'^/2. In Fig. 7(b), we show the extracted 
I{pr,Or) at the angle of 6r — 160° from the right-side 
BRR in Fig. 2. The wave packet extracted from SFA2 
is essentially identical to that from the TDSE except for 
a small shift of the center from 1.18 to 1.21. Here, we 
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FIG. 8: (Color online) Distributions of the momentum of the 
right-side BRR electrons for atomic hydrogen by 5-cycle laser 
pulses with the wavelength of 800 nm at the peak intensities 
of (a) 5.0xlO^^ (b) 1.0x10^'', (c) 2.5xl0", and (d) 5.0 xlO" 
W/cm^ 



see that the radius of BRR is determined by the center 
of the wave packet. 

According to the classical or semiclassical theory the 
maximal energy of the returning electron is 3.17L'p. Ex- 
tending this idea for each optical cycle, we then have 
Er = = 3.17C/p, where Up = A^/A. Consequently, 
we get 



= 1.26 A 



(32) 



For a 5-cycle laser pulse at the peak intensity of 1.0 x 10^^ 
W/cm^ with the wavelength of 800 nm, Ar — 0.93 for the 
right-side BRR, corresponding to pr = 1.17, which agrees 
with the SFA2 calculation. 

The wave packet extracted from ([3T|) should not de- 
pend on 9r if the electron on the BRR comes entirely from 
the backscattering when the scattering angle is large. We 
found this is the case for 6r > 120°, as plotted in Fig. 
7(c). 

We next check the relation ([32)1 for different intensities. 
In Fig. 8 are presented the wave packets extracted from 
the right-side BRR for the photoelectrons generated by 
5-cycle, 800nm laser pulses for peak intensities of 0.5, 1.0, 
2,5 and 5.0 x 10^"' W/cm^. First we note that the wave 
packets derived from the two scattering angles are nearly 
the same. For these pulses, the vector potentials at the 
return time have values Ar = 0.65, 0.93, 1.46 and 2.08, 
respectively. From ([5^ . the radii are calculated to be 
Pr = 0.82, 1.17, 1.84 and 2.62, which compare well with 
the peak positions of 0.86, 1.18, 1.84 and 2.60, from the 
extracted electron wave packets shown in Fig. 8. 



Returning electron wave packet for long pulses 

The analysis of BRR so far has been limited to short 
pulses. For longer pulses, the BRR electrons on each side 
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FIG. 9: Ionization of hydrogen by laser pulses at the peak 
intensity of 1.0 x lO" W/cm^ with the wavelength of 800 
nm and the full durations of 10 and 20 cycles. (a,b): Com- 
parison of the first Born elastic scattering cross sections with 
the SFA2 angular distributions of the BRR electrons; (c,d): 
Distributions oi I{pr,dr) at 6^ = 140° and 180°. 



FIG. 10: Distributions of I{pr, dr) of the BRR electrons, born 
at the time (— 0.5r, — 0.4r) and returning at the Ist-return 
time (0.0, 0.5r), the 2nd-return time (0.5r, l.Or) and the 
3rd-return time (l.Or, 1.5r), at Gr = 170° for ionization of 
hydrogen by a 25-cycle laser pulse at the peak intensity of 
1.0 X 10" W/cm^ with the wavelength of 800 nm. 



can be generated from several born times, each separated 
from the previous one by one full optical cycle. Interfer- 
ence from these coherent electron bursts results in the 
characteristic ATI peaks. 

In Fig. 9, we show the SFA2 results for 10- and 20- 
cycle pulses at the peak intensity of l.Ox 10"'^"' W/cm^ and 
wavelength of 800 nm. For these "long" pulses, accurate 
TDSE results are difficult to obtain. Figs. 9(a) and (b) 
show that the electron yields on the BRR are oscillatory, 
but the envelope is well reproduced by the differential 
elastic scattering cross sections. The oscillations arc at- 
tributed to the electron wave packets, see Figs. 9(c) and 
(d) . Note that the envelopes of the electron wave packets 
in these two figures are similar to the smooth wave packet 
for the 5-cycle pulse shown in Fig. 7(b). The peak posi- 
tions of the envelopes for 10- and 20-cycle pulses are at 
1.16 and 1.17, respectively from the figure, as compared 
to 1.17 from Fig. 7(b). 

For long pulses, electrons which have been released ear- 
lier may return at different times. So far we have con- 
sidered the dominant first return. For later returns it is 
generally believed that the probabilities would be smaller 
since the electron wave packet expected to spread in time. 
Based on the SFA2, we can estimate the wave packets for 
the first and the later returns. For this purpose, we con- 
sider a 25-cycle pulse. We assume that the electrons are 
born in the interval [— 0.5r, — 0.4r]. For the ffist return 
we isolate backscattering occurring within the time in- 
terval [0.0,0.5r]. For the second and the third returns, 
the intervals are chose to be [0.5t, l.Or] and [l.Or, 1.5t], 
respectively, i.e., each is half an optical cycle later from 
the previous return. From the calculated 2D spectra for 
each return, we extracted the electron wave packets, see 
Fig. 10. It is clear that the wave packet from the first re- 
turn has the highest yield and highest momentum. The 



wave packet from second return has lower momentum 
than from the third return, but the peak intensity is 
about twice higher. The wave packets "derived" from this 
model are qualitatively similar to those obtained from the 
classical simulation by Tong et al [30| . but with larger 
width and smaller strength since spreading of the wave 
packet is included in the SFA2 but not in the classical 
simulation of the latter. 



SUMMARY 

In this paper we studied the two-dimensional momen- 
tum spectra of high-energy photoelectrons of a hydrogen 
atom in an intense laser pulse. We focussed on electrons 
that are on the back rescattered ridges (BRR). These 
electrons have been identified initially from solving the 
time-dependent Schrodinger equation (TDSE) and they 
were interpreted as due to the backscattering of the re- 
turning electrons by the target ion. Using the second or- 
der strong field approximation (SFA2) where the rescat- 
tering is accounted for to the first order, we showed that 
BRR electrons also appear in the SFA2 calculation. By 
analyzing the results from SFA2 we have been able to 
identify the time of birth of the electrons and the re- 
turn time where these electrons are backscattered by the 
target ion. From the yields of the electrons along the 
BRR, we further showed that it is possible to extract 
the returning electron wave packet as well as the elastic 
scattering cross sections between the returning electrons 
and the target ion. The electron wave packets extracted 
from SFA2 have been shown to be very close to those ex- 
tracted from solving the TDSE. Since the SFA2 calcula- 
tion is much simpler, this allows us to analyze the return- 
ing electron wave packets conveniently for laser pulses of 
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different intensities and durations. We comment that 
the yield along the BRR calculated using the SFA2 does 
not predict the correct elastic scattering cross section be- 
tween the electron and the target ion since it was based 
on the first order theory. For the present atomic hy- 
drogen target, on the other hand, the elastic scattering 
cross sections calculated from the first-order theory and 
the exact results are identical - a special property of the 
Coulomb potential. While results have been presented 
here only for atomic hydrogen target, similar comparison 
with the same conclusion has been made for other atoms. 
In addition, there is no reason to expect that the same 
conclusion would not hold for molecular targets where 
accurate TDSE solutions at the same level of accuracy 
as in atoms is computationally not feasible. 
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